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摘要 : 非 合作 目标 的 激光 测 距 预报 一 般 是 基于 双 行 根 数 CILE). 外 推出 来 的 ， 往 往 有 较 
大 偏差 ， 对 激光 测 距 的 成 功率 有 较 大 影响 。 结 合 空间 目标 的 轨道 理论 和 实测 时 的 数据 分 析 ， 
预报 的 偏差 主要 由 于 通过 预测 模型 外 推 的 空间 目标 在 运行 轨道 上 的 平 近 点 角 与 实际 平 近 点 
角 存 在 偏差 。 根 据 非 合作 目标 在 望远镜 跟踪 视 场 中 的 脱 靶 量 ， 利 用 相关 算法 可 以 找到 一 个 最 
佳 的 时 间 根 数 偏差 量 修正 空间 目标 的 平 近 点 角 , 经 过 修正 ,空间 目标 的 视 位 置 偏 差 得 到 改善 ， 
距离 偏差 能 够 从 几 百 米 减 小 到 几 十 米 ， 提 升 了 预期 回 波 到 达 时 刻 的 准确 度 ， 可 以 给 单 光 子 探 
测 器 提供 更 高 精度 的 距离 门 控 ， 提 高 测 距 成 功率 。 
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激光 测 距 技 术 是 一 项 综合 技术 , 涵盖 激光 、 电 子 、 天 文 测 量 和 卫星 轨道 计算 等 多 个 学 科 
领域 ,在 实际 工作 中 是 一 个 非常 复杂 的 过 程 ， 受 目标 特性 、 望 远 镜 指向 误差 、 距 离 、 接 收 孔 
径 大 小 等 多 种 因素 的 影响 。 在 众多 影响 因素 中 , 空间 目标 的 轨道 预报 准确 性 和 稳定 性 直接 影 
响 测 距 望 远 镜 发 出 的 激光 脉冲 能 否 持续 到 达 其 表面 , 以 及 接收 端的 单 光子 探测 器 能 否 及 时 探 
测 到 回 波光 子 ， 对 激光 测 距 的 成 功率 有 很 大 影响 。 
卫星 激光 测 距 系统 广泛 采用 了 距离 门 控 技术 ， 即 时 间 滤 波 技 术 ， 有 效 防止 噪声 的 干扰 。 
也 就 是 在 预期 回 波 到 达 时 刻 之 前 一 瞬间 给 单 光子 探测 器 送 开门 脉冲 , 这 样 就 使 得 受到 噪声 干 
扰 的 概率 减 到 最 小 ， 使 得 淹没 在 众多 噪声 中 的 单 光 子 探测 成 为 可 能 凸 。 

非 合 作 目 标的 激光 测 距 预报 ， 即 载 入 跟踪 望远镜 的 轨道 预报 一 般 基 于 双 行 轨道 根 数 
(Two-line elements, TLE)， 经 过 SGP4 (Simplified General Perturbations Version 
4) 轨道 衍生 器 外 推出 来 的 说， 在 对 目标 进行 跟踪 时 往往 有 较 大 的 偏差 ， 主 要 体现 在 望远镜 
视 场 中 的 位 置 偏差 以 及 测 站 与 目标 之 间 的 距离 偏差 , 视 位 置 偏差 过 大 导致 激光 脉冲 不 能 击 中 
目标 ; 测 站 和 目标 的 预报 距离 偏差 过 大 , 使 预期 回 波 到 达 时 刻 与 实际 回 波 到 达 时 刻 出 现 大 的 
偏差 ， 导 臻 单 光子 探测 器 的 距离 门 开启 期 间 探测 不 到 回 波光 子 。 

本 文 提出 的 非 合作 目标 激光 测 距 预 报 实时 修正 方法 的 思路 是 : 先 利 用 双 行 轨 道 根 数 原始 
预报 对 可 见 的 非 合作 目 标 跟踪 几 分 钟 , 获取 其 在 望远镜 跟踪 视 场 中 的 脱 邯 量 数据 信息 , 根据 
这 些 信息 解 算 时 间 根 数 偏差 .并 基于 该 时 间 根 数 偏差 对 轨道 预报 进行 修正 后 , 重新 载 入 望 远 
镜 继 续 跟 踪 目 标 。 激 光 测 距 卫 星 (例如 ajisai、beaconc、starlette、stella 等 ) 具有 双 行 
轨道 根 数 和 CPF(Consolidated Prediction Format) 两 种 预报 ， 文 [5] 将 中 低 轨 卫星 的 CPF 
星 历 与 利用 全 球 激光 跟踪 站 的 观测 数据 定 轨 出 的 精密 轨道 进行 了 比较 , 分 析 CPF 星 历 精度 在 
米 级 以 内 ， 文 [6] 分 析 了 用 SGP4 模型 处 理 近 地 目标 时 ， 定 轨 精 度 在 百 米 量 级 。 为 了 排除 望 远 
镜 指向 误差 等 测 距 系 统 误差 的 影响 , 在 理论 上 证 明 该 方法 的 有 效 性 , 本 文 把 中 低 轨道 高 度 的 
激光 测 距 卫 星 视 为 非 合作 目标 进行 研究 分 析 , 把 某 一 时 间 段 可 见 的 激光 星 的 CPF 预报 作为 参 
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考 值 ， 即 模拟 观测 值 。 目 标的 脱 靶 量 由 同一 时 刻 双 行 轨道 根 数 预报 和 CPF 预报 中 的 方位 角 
和 高 度 角 侦 差 量 代替 ， 然 后 进行 时 间 根 数 侦 差 的 解 算 和 轨道 预报 的 修正 。 


1 激光 测 距 预报 的 生成 
1. 1 时 间 系 统 的 转换 


协调 世界 时 〈UTC) 是 一 种 修正 过 的 原子 时 ， 空 间 目 标的 状态 以 及 观测 数据 都 使 用 UTC 
作为 输入 、 输 出 的 时 间 参 数 。 格 林 尼 治平 恒星 时 (Greenwich sidereal time ,GST) 96 是 
坐标 系 转换 的 一 个 关键 参数 ， 故 要 由 协调 世界 时 计算 出 9e， 详 见 文 [7]。 


1.2 TEME 坐标 系 与 站 心地 平 坐标 系 之 间 的 转换 


利用 预测 模型 外 推出 来 的 空间 目标 位 置信 息 是 在 真 赤道 平 春分 点 坐标 系 (True Equator 
Mean Equinox Frame, TEMEO 中 表达 的 ， 而 激光 测 距 望远镜 使 用 的 轨道 预报 则 是 在 站 心 
地 平 坐标 系 中 表达 的 ,所 以 要 实现 由 真 赤 道 平 春分 点 坐标 系 向 站 心地 乎 坐标 系 的 转换 ， 如 图 
1. 


图 1 TEME 坐标 系 向 站 心地 平 坐标 系 转换 示意 区 
Fig.1 Schematic diagram of the conversion from the TEME system to the station local system. 
cosg sing 0 
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RIR Hat ECASGÉCE- 484) SUA RR MERE: Rz(96) 为 转换 矩阵 ，06 为 格林 尼 
治平 恒星 时 ，7 为 空间 目标 在 地 固 坐标 系 中 的 位 置 矢量 ， 广 为 空间 目标 在 站 心地 平 坐标 系 中 
的 位 置 矢 量 ; Fs** 为 激光 测 距 望远镜 在 地 固 坐 标 系 中 的 位 置 矢量 ;，(MLT) 为 转换 矩阵 : 


(MLT) = R, E E $) R, E + i) ; (2) 
1 0 0 
其 中 ， R.($- 9)» ( sing cose) (3) 
0 —cosp sing 


—sinA cosà 0 
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e (2) = (a —sinA 0 


$ 和 4 为 激光 测 距 望远镜 所 在 位 置 的 纬度 和 经 度 ，R (Z o MR, (Z + ) 均 为 转换 矩阵 。 


1.3 光 行 差 改正 


当 对 中 低 轨 空间 目标 进行 激光 测 距 时 , 测 站 和 目标 之 间 往 往 具有 几 千 干 米 的 距离 ,在 激 
光 脉 冲 的 上 行 和 下 行 期 间 , 空间 目标 相对 测 站 发 生 了 较 大 变化 , 因此 对 于 激光 测 距 一 般 用 发 
射 时 刻 作为 时 标 生 成 跟踪 文件 ， 即 望远镜 的 实时 指向 有 个 前 指 量 ,而 针对 本 文 的 应 用 ， 由 于 
下 行 光 行 差 的 原因 , 目标 可 见 对 应 的 预报 应 该 滞后 , 即将 激光 接收 方向 作为 模拟 观测 值 。 CPF 
预报 中 ， 前 指 量 的 算法 如 下 : 
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对 于 某 一 时 刻 T， 空 间 目标 和 测 站 在 地 固 坐 标 系 中 的 位 置 矢 量 为 (X,Y,2Z) 和 ?sto (X, 
加 ，20)， 则 测 站 和 空间 目标 之 间 的 相对 距离 D = [P 一 产 t?|。 在 不 考虑 光 行 差 的 情况 下 ， 激 
光 脉 冲 发 射 方向 即 T 时 刻 空间 目标 的 预报 位 置 ， 激 光 脉 冲 在 上 行 光路 传播 的 时 间 Tojt = 
D/c ，c 为 光速 ， 空 间 目 标 在 站 心地 平 坐标 系 中 的 位 置 矢量 (px, py,pz)。 在 考虑 光 行 差 的 
情况 下 ， 激 光 脉 冲 发 射 方向 为 Ts =T + Tout 时 刻 空间 目标 的 预报 位 置 。 空 间 目 标 在 站 心地 
平 坐标 系 中 的 位 置 矢量 (psx Pey pgz)。 此 处 需要 特别 注意 的 是 地 固 坐 标 系 是 一 个 旋转 坐 
标 系 ， 计 算 中 需要 进行 相应 改正 。 

这 样 就 得 到 上 行 光路 的 光 行 差 引起 的 激光 发 射 前 指 量 : 
dA=A—Ap=tant (- 2 — tan“! (- m) ; (5) 


x 


dE = E — Ep = tant L — tant 2 , (6) 
E 55 px?* py? e DBx^ * pBy? 
XIF OPBATOGERBUJC6TT 28 5 e PH OG Belit Jer s RIF ZR TRI J LER AAA READS R23 
dA, dE. 


—— 


2 时 间 根 数 偏差 的 解 算 以 及 预报 的 修正 


2.1 理论 表述 


环绕 地 球 做 椭圆 运动 的 空间 目标 的 轨道 示意 图 如 图 2。 


K 
Earth's north polar axis 
e 


Node line 


图 2 环绕 地 球 运 行 的 空间 目标 的 轨道 示意 
Fig2 Schematic diagram representing the orbit of earth-circling spatial target. 


图 2 中 , 阴影 部 分 为 地 心 天 球 坐 标 系 O-XYZ 的 基本 平面 Q^ rl) , 接近 J2000.0 平 赤道 
面 ，X 轴 指向 接近 J2000. 0 平 春 分 点 ，Z 轴 指向 地 球 的 北极 。 7 为 空间 目标 的 轨道 半径 ，h 为 
空间 目标 轨道 平面 的 法 向 量 ， 尺 为 通过 轨道 升 交点 的 矢量 ，8 为 通过 轨道 上 近地点 的 矢量 。 
为 从 了 轴 方 向 起 量 的 轨道 升 交点 方向 六 的 经 度 ;7 为 轨道 面 与 基本 平面 的 夹 角 ; oo 为 近地点 
幅 角 ;6 为 极 坐标 变量 。 空 间 
_ all-e’ (T) 


1-ecos(0—o) ` 

E CO 式 中 ，o 为 半 长 径 ，e 为 偏心 率 。 近 地 点 幅 角 w 和 极 坐 标 变量 9 都 是 从 升 交 点 方向 

起 量 的 ， 当 有 摄 动 时 ， 椭 圆 轨 道 随时 间 变 化 ， 升 交点 方 同 也 在 变化 ，@ 应 从 这 变化 的 升 交点 
方向 起 量 ， 而 极 坐标 变量 6 仍 从 一 个 定义 的 不 变 方向 起 量 。 


M -n(t- 1) ; (8) 
从 近地点 开始 起 量 ， 其 中 ， 为 平 运动 角速度 。 空 间 目 标的 平 近 点 角 是 一 个 与 时 间 相关 的 参 
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数 ， 能 够 确定 对 应 时 刻 目 标 在 轨道 面 内 椭圆 轨道 上 的 位 置 。 当 上 = Tt 时 ，IJE0， 相 应 地 ， 空 间 
目标 在 近地点 处 ，r 达 到 极 小 值 ，Tt 就 是 空间 目标 过 近地点 的 时 刻 。 

空间 目标 在 轨 运 行 时 会 受到 多 种 摄 动力 的 影响 ， 包 括 地 球 的 非 球形 引力 、 日 月 引力 、 大 
气 阻力 ， 太 阳 辐 射 压力 等 ， 对 于 处 于 中 低 轨道 的 空间 目标 ， 受 到 大 气 阻力 的 影响 比较 大 。 双 
行 轨道 根 数 考虑 了 地 球 非 球形 引力 、 日 月 引力 的 长 期 和 周期 摄 动 影响 ,以 及 大 气 阻 力 模型 产 
生 的 引力 共振 和 轨道 衰退 。 双 行 轨道 根 数 是 “平均 ” 根 数 , 它 用 特定 的 方法 去 掉 了 周期 扰动 
项 , 预测 模型 SGP4 用 同样 的 方法 重 构 了 周期 扰动 项 中。 在 观测 中 ， 空 间 目 标 往往 在 沿 迹 方向 
与 望远镜 的 指向 中 心 有 较 大 偏差 ， 提 前 或 者 滞后 ， 主 要 原因 是 SGP4 模型 采用 的 大 气 模 型 简 
化 为 B 一 个 阻力 参数 ， 沿 迹 方 向 的 偏差 量 可 以 简化 归 因 于 时 间 根 数 t 偏 差 。 

如 果 在 跟踪 目标 时 ， 适当 调整 t 的 值 ， 可 以 减 小 空间 目标 与 望远镜 指向 中 心 的 偏 移 量 , 修 
正 空 间 目 标 与 测 站 的 距离 量 ， 并 提高 跟踪 的 稳定 性 。 

时 间 根 数 t 偏 差 的 解 算 算法 如 下 ， 记 空间 目标 在 站 心地 平 坐标 系 中 的 方位 角 为 4， 高 度 
角 为 6， 离 测 站 距离 值 为 s， 实 测 的 序列 记 为 (ti, Ap Eps) i= 1,2,3...n， 在 计算 预报 时 适当 
调整 rz 的 值 ， 即 令 z + Ar 取代 r，Az 是 调整 量 ， 亦 即 空间 目标 在 真 赤 道 平 春分 点 坐标 系 中 相应 
IF] 22] ER BER OX Yo Zi ti + ADAR R(X; Yi, Zo ti): 

F(t) = R/(05())R(t AT), BCE) = (MLT) — 7*9]; 
由 此 可 得 预报 序列 (ti, A'i, E'u s'u Av) i= 1,23 n o WAA; = A'; — Ap AE; = E'; — Ep At 


I . ; à ,2 
RMS* [Zea A (9) 


为 跟踪 时 间 段 上 的 极 小 值 ， 此 时 Ar 即 为 预报 的 时 间 根 数 偏差 。 上 述 问 题 是 一 个 复杂 的 非 线 
性 问题 ， 直接 求解 比较 困难 ， 考 虑 到 时 间 根 数 偏差 精确 到 毫秒 即 可 ， 因 此 一 个 简单 直接 的 方 
式 就 是 在 一 定 范围 内 进行 搜索 。 

在 空间 目标 的 指向 信息 得 到 修正 时 , 距离 信息 也 得 到 了 相应 的 修正 , 进而 能 够 提高 单 光 
子 探测 器 的 距离 门 控 的 精度 ， 增 大 接收 回 波 光子 的 概率 。 


2. 2 实例 分 析 

ajisai、beaconc、starlette、stella 等 激光 测 距 卫星 位 于 近 地 轨 道 ， 与 大 部 分 的 非 合 作 
目标 的 轨道 高 度 相仿 , 空间 环境 类 似 , 能 够 很 好 地 替代 仅 有 双 行 轨道 根 数 预报 的 非 合作 目标 
进行 研究 分 析 。 

本 文 以 Stella. 激光 测 距 卫星 的 过 境 情 况 进 行 分 析 。Stella 卫星 的 轨道 高 度 为 
792.3X803.1km， 处 于 非 合作 目标 的 密集 区 。 当 对 空间 目标 进行 激光 测 距 时 ， 测 站 和 目标 之 
WERA ETTER, 在 激光 脉冲 的 上 行 和 下 行 期 间 , 空间 目标 相对 测 站 发 生 了 较 大 
变化 , 指向 滞后 量 如 图 3, 横 轴 为 日 积 秒 , 纵 轴 分 别 为 方位 角 偏 差 和 高 度 角 偏差 。 可 以 看 出 ， 
滞后 量 波 动 幅 度 比较 大 ， 最 高 时 可 达 十 几 角 秒 ， 是 解 算 时 间 根 数 偏差 时 不 可 忽视 的 因素 。 
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图 3 激光 星 Stella 某 一 过 境 弧 段 中 望远镜 接收 回 波光 子 应 有 的 指向 滞后 量 。 
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Fig3 The lag amount telescope should remedy when receiving the echo photon in one passing segment of Stella 
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图 4 激光 测 距 卫星 stella 分 别 以 过 境 弧 段 的 前 1. 2. 3 分 钟 和 整个 弧 段 的 跟踪 数据 搜索 出 的 时 间 根 数 偏 差 结 


100 150 200 


Fig.4 The time element deviations of stella calculated by using the tracking data of previous one. two. three 


minutes and entireness of the passing segment. 


图 4 显示 了 在 stella 的 某 次 过 境 观测 弧 段 ， 分 别 用 前 1min, 2min, 3min 以 及 全 弧 段 跟踪 
数据 进行 时 间 根 数 搜索 的 结果 ， 图 中 横 轴 为 时 间 根 数 偏差 Ar ， 纵 轴 为 对 应 的 RMS*， 时 间 根 
论 用 前 1min. 2min, 3min 还 是 全 弧 段 跟踪 数据 进行 
搜寻 都 有 一 个 极 小 值 ， 对 应 时 间 根 数 偏差 分 别 为 26ms、27ms、29ms、32ms， 并 无 显著 差 


数 搜寻 范围 为 -100,200]。 由 图 可 见 ,， 无 


异 。 对 于 其 它 的 激光 测 距 卫 星 , 根据 不 同时 长 的 跟踪 弧 段 , 搜索 的 时 间 根 数 偏差 结果 如 表 1. 
表 1 若干 激光 测 距 卫 星 根 据 不 同时 长 的 跟踪 数据 解 算出 的 时 间 根 数 偏差 
Table 1 The time element deviations of some Laser Ranging Stallites by using 
different duration of the tracking data. 


激光 星 轨道 高 度 (km) 过 境 时 长 
ajisai 1479X1497 12min20s 
beaconc 928.6X1305 8min30s 
starlette 805.5X1107 6min9s 


前 一 分 钟 (ms) 


35 


30 


36 


前 两 分 钟 (ms)  — 前 三 分 钟 (ms) 全 弧 段 (ms) 
35 35 34 
28 26 23 
38 40 43 


通过 表 1 可 以 看 出 ,上述 3 个 激光 测 距 卫星 在 过 境 可 见 的 时 间 段 内 , 用 前 1min、2min、 


3min 以 及 整个 弧 段 解 算 的 时 间 根 数 偏差 无 显著 差异 。 根 据 搜索 时 间 根 数 偏 差 的 思想 可 知 ， 


利用 的 跟踪 数据 段 越 长 , 解 算 的 结果 越 精确 , 不 


| 用 全 弧 段 跟踪 数据 解 算 的 结果 在 理论 上 是 最 


精确 的 。 考 虑 到 中 低 轨 目标 过 境 时间 短 等 问题 ， 选 择 前 1min 数据 解 算 的 时 间 根 数 偏差 适用 


于 整个 弧 段 是 比较 理想 的 。 
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图 5 stella 的 轨道 预报 分 别 以 Oms 〈 未 修正 )、26ms、27ms、29ms、 
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32ms 为 时 间 根 数 偏差 修正 后 与 参考 轨道 在 高 度 角 Ca) 和 方 


Fig5. The deviations of stella’ orbit prediction corrected by Oms, 26ms, 27ms 29ms 32ms compared with the reference 


orbit in azimuth angle and elevation angle. 


22824 (20180311 20:37:30.000-20:44:00.000,792.3X803.1km) 


T T 
pred( Oms)-cpf 
pred(26ms)-cpf 


AE/ 


-30 L 1 1 1 L 1 L 
0 20 40 60 80 100 120 140 160 


AA*cosE/ 


图 6 stella 的 轨道 预报 分 别 以 Oms 〈 未 修正 )、26ms 为 时 间 根 数 偏差 修正 后 与 参考 轨道 的 偏差 量 。 


Fig6. The deviations of stella’ orbit prediction corrected by Oms, 26ms compared with the reference orbit. 
激光 卫星 stella 的 轨道 预报 经 过 相应 的 修正 后 ， 视 位 置 偏差 明显 减 小 ， 且 波动 幅度 也 乱 

小 。 图 5(a) 中 ， 实 线 表 示 stella 的 TLE 原始 预报 中 的 方位 角 与 CPF 参考 值 的 残 差 随 观 测 时 间 
A 的 变化 ， 下 方 的 虚线 分 别 表示 TLE 预报 经 过 26ms、27ms、29ms、32ms 修正 的 方位 角 与 CPF 
© 参考 值 的 残 差 随 观测 时 间 的 变化 。 类 似 地 ， 图 5(b) 中 ， 实 线 表示 stella 的 双 行 轨道 根 数 原始 
预报 中 的 高 度 角 与 CPF 参考 值 的 残 差 随 观测 时 间 的 变化 , 中间 的 虚线 分 别 表示 双 行 轨道 根 数 
预报 经 过 26ms、27ms、29ms、32ms 修正 的 高 度 角 与 CPF 参考 值 的 残 差 随 观测 时 间 的 变化 。 
CA 图 6 是 以 AA * cosE 为 横 坐 标 ，AE 为 纵 坐 标 画 出 了 整个 跟踪 弧 段 以 Ooms (未 修正 ) 和 26ms 
T (前 一 分 钟 ) 分 别 修正 后 的 视 位 置 偏 差 量 的 变化 。 综 合 以 上 可 以 看 出 ，TLE 预报 经 过 时 间 根 
C2 数 偏差 改正 后 ,， 视 位 置 偏差 量 明 显 减 小 ， 且 在 某 一 固定 值 附近 小 范围 波动 , 波动 范围 在 几 角 

秒 之 内 。 

表 2 若干 激光 测 距 卫星 修正 后 的 轨道 预报 在 方位 角 、 高 度 角 的 偏差 均值 以 及 相应 的 RMS*。 


Table 2 The mean value of the deviations of some Laser Ranging Stallites' corrected 


orbit predictions in azimuth angle and elevation angle and corresponding RMS*. 


激光 星 AA, cos Ep AE, AA, cos E, AE, AA; cos E; AE, RMS RMS RMS} 

ajisai -14.7301 15.1602 -20.1686 4.0037 -20.0124 4.3220 19.2066 13.5802 13.5826 
beaconc -13.4514 4.0645 -9.3972 -8.8873 -10.3501 -5.8718 12.1740 8.7697 4.1602 
starlette -69.3258 15.4581 -62.3066 -0.3636 -60.9365 -3.4399 89.6907 84.2076 87.8905 
stella 13.2867 18.5682 4.5302 5.7109 2.5032 2.7460 13.6939 1.7213 2.4036 


表 2 中 展示 了 处 于 近 地 轨 道 的 4 颗 激 光 卫 星 的 双 行 轨道 根 数 预 报 在 未 修正 及 基于 前 一 分 
钟 ， 整个 弧 段 解 算出 的 时 间 根 数 偏 差分 别 修正 后 与 CPF 参考 轨道 在 方位 角 、 高 度 角 的 偏差 均 
值 以 及 相应 的 位 置 偏差 RMS*。AAo cos E60，AEo 是 双 行 轨道 根 数 原始 预报 与 参考 轨道 在 方 
位 角 、 高 度 角 在 整个 过 境 弧 段 的 偏差 均值 。AA1 cosE AE.: A42 cosE2，AE2 分 别 是 基于 
前 一 分 钟 、 整 个 弧 段 解 算 的 时 间 根 数 偏差 进行 修正 后 的 轨道 预报 与 CPF 参考 轨道 在 方位 角 、 
高 度 角 的 偏差 均值 。 可 以 看 出 后 两 者 的 位 置 偏差 RMS* 比 双 行 轨道 根 数 原始 预报 的 要 小 ， 即 
更 接近 实测 运行 轨迹 。 后 两 者 在 方位 角 、 高 度 角 的 偏差 均值 也 很 接近 ,更 进一步 说 明 可 以 用 
前 一 分 钟 的 跟踪 数据 解 算 出 的 时 间 根 数 偏差 代 蔡 整个 弧 段 的 跟踪 数据 解 算 的 时 间 根 数 偏差 
进行 轨道 预报 修正 。 
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图 6 激光 星 stella 的 轨道 预报 以 不 同 的 时 间 根 数 偏差 修正 后 ， 距 离 预报 值 与 参考 值 的 偏差 量 。 


Fig6. The distance deviations of stella' orbit prediction corrected by different time elememt bias compared with 


reference values 
ps 从 图 6 可 以 看 出 ， 激 光 星 stella 的 距离 预报 分 别 以 26ms、27ms、29ms、32ms 为 时 间 根 
C 数 偏差 修正 后 ， 距 离 偏差 由 最 大 的 二 百 多 米 减 小 到 几 十 米 ， 并 且 变 化 范围 更 稳定 。 
表 3 若干 激光 测 距 卫星 修正 后 的 轨道 预报 在 预报 距离 方面 的 偏差 均值 以 及 相应 的 RMS。 


Table 3 The mean value of the deviations of some Laser Ranging Stallites' corrected 


orbit prediction in range and corresponding RMS 


激光 星 Do n D, RMS) RMS, RMS, 
Ajisai -215.89 -49.73 -54.52 78.80 9.68 11.60 
Beaconc -154.19 8.72 -29.21 46.16 9.11 4.35 
Starlette -214.18 -4.56 36.19 56.31 14.84 7.67 
stella -180.97 -37.51 -4.99 62.42 18.58 8.67 


表 3 展示 了 处 于 近 地 轨 道 的 4 颗 激 光 卫 星 的 轨道 预报 在 未 修正 及 基于 前 一 分 钟 ， 整 个 弧 

段 解 算出 的 时 间 根 数 偏差 分 别 修 正 后 相 较 于 CPF 参考 轨道 在 预报 距离 方面 的 偏差 均值 Do、 

D1、D, 以 及 相应 的 RMS， 即 RMSo、RMS1、RMS,。 可 以 看 出 ， 预 报 距 离 偏 差 均 值 均 由 几 百 

米 减 小 到 了 几 十 米 ， 甚 至 几米 ;， 均 方 根 值 RMS 也 减 小 了 一 个 量 级 ， 说 明 偏差 稳定 性 显著 增 

(S) 强 。 上 述 结果 说 明基 于 时 间 根 数 偏差 修正 双 行 轨道 根 数 预报 能 够 提高 单 光 子 探测 器 的 距离 门 
控 的 精度 。 


结论 


本 文 介绍 了 一 种 非 合 作 目 标 激光 测 距 预 报 实时 修正 的 方法 。 在 搜索 时 间 根 数 偏 差 的 过 程 
中 , 以 RMS* 取 得 极 小 值 时 的 Ar 值 为 时 间 根 数 偏差 ,用 于 空间 目标 的 双 行 轨道 根 数 预报 修正 ; 
用 不 同时 长 的 跟踪 数据 搜寻 出 的 时 间 根 数 偏差 无 显著 差异 ,用 前 一 分 钟 数 据 解 算 的 时 间 根 数 
偏差 可 以 适用 于 整个 弧 段 ; Se 空间 目标 的 实测 运行 轨迹 且 预 
报 距离 偏差 值 均 由 几 百 米 减 小 到 了 几 十 米 , 甚至 几米 ,能 够 提高 单 光 子 探测 器 的 距离 门 控 的 
精度 。 

本 文选 取 的 若干 个 空间 目标 的 双 行 轨道 根 数 预报 经 过 此 方法 的 修正 后 ， 视 位 置 值 和 预报 
距离 值 的 偏差 得 到 了 明显 的 改善 。 
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Research on Real-time correction method of laser 


ranging prediction of noncooperative target 


Zhang Xunfang!?, Zhao Xue!?, Li Rongwang!?, Li Zhulian'? 
(1.Yunnan Observatories,Chinese Academy of Sciences, Kunming 650011,China; 
2.University of Chinese Academy of Sciences,Beijing 100049,China; 
3.Key Laboratory of Space Object and Debris Observation, Chinese Academy of 
Sciences,Nanjing,210008) 

Abstract: The laser ranging prediction of noncooperative target is extrapolated from TLE 
commonly, existing large deviation, having great influence on the mission success rate of laser 
ranging. Based on the orbital theory and the measured data of space target, the deviation of the 
forecast is mainly due to the deviation of mean anomaly from the prediction model to the real 
mean anomaly of the spatial target in the orbit. According to the target's miss distance in the field 
of view of tracking telescope, the relative algorithm can be used to find the optimal time element 
deviation to correct the mean anomaly of the spatial target. After modification, the apparent 
position deviation of the spatial target is decreased, the distance deviation can be reduced from 
several hundred meters to tens of meters making the accuracy of expected echo arrival time 
better,providing the more precise distance gate control for the single-photon detectorimproving 
the mission success rate of laser ranging. 

Key words: real-time correction; time element deviation; noncooperative target; orbit prediction 
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